Priedmann cosmology with a generalized equation of state and bulk viscosity 
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The universe media is considered as a non-perfect fluid with bulk viscosity and described by a 
more general equation of state. We assume the bulk viscosity is a linear combination of the two 
terms: one is constant, and the other is proportional to the scalar expansion 8 = 3a/a. The equation 
of state is described as p — (7 — l)p + po, where po is a parameter. This model can be used to 
explain the dark energy dominated universe. Different choices of the parameters may lead to three 
kinds of fates of the cosmological evolution: no future singularity, big rip, or Type 111 singularity of 
Ref. [S. Nojiri, S.D. Odintsov, and S. Tsujikawa, Phys. Rev. D 71, 063004 (2005)]. 



PACS numbers: 98.80.Cq, 98.80.-k 
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I. INTRODUCTION 

The cosmological observations indicate that the ex- 
pansion of our universe accelerates 0]. Recently lots of 
work on extended gravity 0]such as modifying equation 
of state or by introducing the so called dark energy is 
to explain the cosmic acceleration expansion observed. 
To overcome the drawback of hydrodynamical instabil- 
ity, a linear equation of state of a more general form, 
p = a(p — po) is proposed 3], and this form incorporated 
into cosmological model can describe the hydrodynami- 
cally stable dark energy behaviors. 

The observations also indicate that the universe media 
is not a perfect fluid and the viscosity is concerned in 
the evolution of the universe 0, IS • ^ n the standard 
cosmological model, if the equation of state parameter us 
is less than -1, the universe shows the future finite singu- 
larity called Big Rip HQ. Several ideas are proposed to 
prevent the big rip singularity, like by introducing quan- 
tum effects terms in the action. 

In this paper, we show that the Friedmann equations 
can be solved with both a more general equation of state 
and bulk viscosity detailed as follows. The equation of 
state is 



P= (7- l)p+Po, 



(1) 



where po and 7 are two parameters. The bulk viscosity 
is expressed as 



C - Co + Ci 



(2) 



where Co and £1 are two constants conventionally. The 
uj = p/p is constrained as —1.38 < u < —0.82 [l(| by 
present observation data, so the inequality in our case 
should be 



4.38 < 7- 1 + — < -0.82. 
P 



(3) 
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The parameter pq can be positive (attractive force) or 
negative (repulsive force), and conventionally £0 and £i 
are regarded as positive. To choose the parameters prop- 
erly, it can prevent the Big Rip problem or some kind of 
singularity for the cosmology model, like in the phantom 
energy phase, as shown below. 

This paper is organized as follows. In Sec. II we de- 
scribe our model and give out the exact solution. In Sec. 

III we consider some special cases of the solution. In Sec. 

IV we discuss the acceleration phase and the future sin- 
gularities in this model, and in the last section (Sec. V) 
we summarize our conclusions. 



II. MODEL AND CALCULATIONS 

We consider the Friedamnn-Roberson- Walker metric 
in the flat space geometry (k=0) 

ds 2 = -dt 2 + a(t) 2 (dr 2 + r 2 dVL 2 ), (4) 

and assume that the cosmic fluid possesses a bulk viscos- 
ity C- The energy-momentum tensor can be written as 



T llv = pU vt U v + {p-^9)H lw , 



(5) 



where in comoving coordinates l7 m = (1,0), 
3(i/ a, and = Qfiu + U[JJ V [ll(. By defining the effec- 
tive pressure as p = p — Q6 and from the Einstein equation 
Rfj,v ^ \g^vR — S^GT^, we obtain the Friedmann equa- 
tions 



8ttG 
4ttG 



{p + 3p). 



(6a) 
(6b) 



The conservation equation for energy, T.°", yields 

p+(p + p)8 = Q. (7) 

Using the equation of state to eliminate p and p, we 
obtain the equation which determines the scale factor 
a(t) 

- = -^-^^ + l^GCo--^G Po , (8) 
a 2 gt a 
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where the effective equation of state parameter is shifted 
from the original one as 

7 = 7-8ttGCi. (9) 

So we can see that the equivalent effect of the second term 
in £ is to change the parameter 7 to 7 in the equation of 
state. As shown in Ref. [5J, the barrier u> = — 1 between 
the quintessence region (u> > —1) and the phantom region 
{uj < —1) can be crossed, as a consequence of the bulk 
viscosity available. 

Since the dimension of the two terms 127rG(o and 
—AirGpo is [time] -1 and [time] -2 , respectively, we define 

12^GCo = ^ (10) 
-AirGpo = i (11) 

I 



then Eq. (jSJ becomes 



a _ 37 - 2 a 2 la 1 
a ~ 2 ^2 + jT~ + f2- 



Here Ti and T 2 are criteria to determine whether 
we should concern the £ and pq. If Ti >> 
t, thecosmictimescale, the effect of £ can be neglected, 
and if Ti >> t, the effect of po can be neglected likewise. 

Concerning the initial conditions of a(t ) = a and 
0(*o) = ^0; if 7 0, the solution can be obtained as 



t-*o 



1 



T 

1®qT + — 1 oq> 



t-tr 



1 

3\ 
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And we obtain directly 



1 (1 + j6 T + ^)exp(^) - (1 - 7 O T + - 



a 3 7 (i + 7 o T-^)exp(^) + (l-70 o r+^) 

I 



(13) 
(14) 



Here we define 



T = 



Ti 



(15) 



We can see that when T2 — ► 00, T = Ti; when Tl — ► 00, 
T = r 2 /V67. 



III. 7 = AND SPECIAL CASES 

For 7 = 0, we should use the mathematical 
L'Hospital's rule to calculate the limit of Eq. ijT3|l rigously 
and note 



lim < ^-— 

7-»o cry 



3T? 

J 2 



(16) 



The limit of solution a(t) when 7 — > is 

a(t) = aoeM^OoT, + ^ (e^)^ ^ 



Ti(t-t ). 



T 2 



Directly solving Eq. (fT2"|) 



1 a 1 

Ti a T 2 



(17) 



(18) 



gives the same result. So Eq. (|13|) is consistent for 7 
crossing zero. The T 2 — > 00 limit of Eq. i|17|) is 



a(t) = aoexp 
and the T\ — > 00 limit of Eq. (|T7|) is 

a(i) = aoexp 



-V {t - t ) H 



(19) 



(20) 



These two special limits are also consistent with directly 
solving Eq. by checking. 

Let us discuss two special cases in the following. When 
the constant term in the equation of state is not con- 
cerned, i.e. T 2 — > 00, 

2/37 



a(t) = a 



1 



(21) 



and when the constant term in the bulk viscosity is not 
concerned, i.e. Ti — > 00, 



a(t) = a { cosht -7£fT- + jOoTsinh ^ 



t-t 



2/37 



When T — ► 00, the two cases become 



a(t) = a 



2/37 



(22) 



(23) 
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For 7 — > 0, the limit case is 

o(t) = aoe 9 ^-^/ 3 , 



(24) 



which corresponds to the de Sitter universe with acceler- 
ating cosmic expansion. 

Additional notions: Eqs. ijfiajl and (|6b|) can be rewrit- 
ten as 



H z = 



8ttG 



H = -4ttG(P + p). 



(25a) 
(25b) 



From these equations, the relation among viscosity, the 
scalar factor a and Hubble parameter H is 

aH^- = -^H 2 + 12ttG(H -6ttG Po . (26) 
da 2 

which reflects the viscosity functions for dark energy and 
matter dominated universe evolution. 



IV. ACCELERATION AND BIG RIP 



If the universe accelerates, then mathematically 



> 0. 



(27) 



From Eq. |(HJ), we can qualitatively see that the bulk vis- 
cosity and a negative density po can cause the universe to 
accelerate. Since the expression of a /a is too complicated 
in this situation, now we only discuss a special case, with 
Po = 0. Here d/a > yields 

= 7-1 > H e (*-*°)/ T - 1 + -L+SttGCi. (28) 

6 U J- 

As we know, if the bulk viscosity is zero as in the standard 
FriedammnRobertson- Walker cosmology model, an ac- 
celerating expansion universe corresponds to lo < —1/3. 
Inequality l|28|l tells us that if the bulk viscosity is large 
enough, the universe expansion can accelerate even if 
w > -1/3. 

According to [9J , the future singularities can be classi- 
fied in the following way: 

• Type I ("Big Rip"): For t -> t s , a -> oo, p -> oo 
and \p\ — > oo 



• Type II ("sudden"): For t -> t s , a 
and \p\ — ► oo 



a s , P 



• Type III: For t — *■ t s , a — > a s , p — >■ oo and \p\ 



oo 



• Type IV: For t — > t s , a — > a s , p — > 0, \p\ — > and 
higher derivatives of H diverge. 

In this paper, p — > oo means p — > oo (we assume 7 / 1 
generally). In the following we show that different choices 
of the parameters may lead to three fates of the universe 
evolution: no future singularity, big rip, or the Type III 
singularity. 



7 < 



From Eq. I|6ajl. we see 



1- a 
a 



(29) 



If the denominator of Eq. (|14fl is zero, 



1 + - Q exp (^^) +1- 7 O T+^ = O. 

(30) 

then a — > 00, p — > 00, so the big rip occurs. The solution 
for t is 



ts = Tin - 



i- 7 fl r + ^' 

l + ^ T-£ 



(31) 



If we want to prevent the big rip, there should be no real 
solution for t > to, so 



1 + WT 



< 1. 



The inequality is equivalent to 

1 + 7 #oT - ^ > 0. 



(32) 



(33) 



The above inequality can be satisfied in some conditions 
for the phantom energy, so the big rip will not occur. 
Furthermore, we can see that even if the dark energy 
is in the quintessence region, there also can be future 
singularity. For example, if 7 > and pq < 0, it is 
possible that inequality is not satisfied, so the future 
singularity may occur, which will be discussed in the next 
subsection below. 

The more explicit form of inequality (|33|l is 



1 



76*0 Ti 



1 



y/1 + 6 7 (Ti/T 2 ) 2 y/1 + 6 7 (Ti/T 2 ) : 



> 0. (34) 



From this inequality, we obtain that 

(i) po < 0, i.e. T| > 0: inequality ffify is always 
unsatisfied, so there will be a big rip at time t s . 

(ii) po > 0, i.e. T| < 0: inequality ijSlfl) is not always 
unsatisfied. If it is unsatisfied, there will be a big rip at 
time t s ; if it is satisfied, there is no future singularity. 



B. 7 > 

If the denominator of inequality l|14|) is zero, then a — > 
0, p — > 00, so the Type III singularity occurs. Following 
the same steps as before, we obtain that 

(i) po < 0, i.e. T| > 0: inequality l|3^|) is always 
satisfied, so there is no future singularity. 

(ii) po > 0, i.e. T| < 0: inequality ffify is not always 
satisfied. If it is satisfied, there is no future singularity; 
if it is unsatisfied, there will be Type III singularity at 
time i=. 
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FIG. 1: pa = 0, 7 = 1, the dash line corresponds to Ti = 
5.1 x 10 28 s, and another curve corresponds to Ti = 5.1 x 10 17 s. 



FIG. 3: £ = 0, 7 = 1, the dash line corresponds to T2 — 5.1 x 
10 16 is, while another curve corresponds to T2 = 5.1 x 10 17 s 
case. 
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FIG. 2: £ = 0, po = 0, the dash line corresponds to 7 = 0.181, 
another curve corresponds to 7 = 0.18 case. 
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FIG. 4: C = 0, T 2 = 5.1 x 10 16 is, 7 = 0.1. 



CONCLUSION 



7 = 



For the case 7 = 0, 



- = \e^-^ + S (e<* 

a 3 i, V 



-to)/T! _ 



(35) 



So there is no future singularity in this case. 

To illustrate the parameters in the general solution Eq 
(|13fl more clearly, we draw some graphics in Fig. 1-4. The 
initial condition is 11] to = 1000s, O = 1.5 x 10 _3 s _1 . 
At this time, the bulk viscosity £ = 7.0 x 10~ 3 g/cm ■ s, 



the corresponding T\ 
for simplicity. 



5.1 x 10 8 s. We assume C,\ 



In conclusion, we have solved the Friedmann equations 
with both a more general equation of state and bulk vis- 
cosity, and discussed the acceleration expansion of the 
universe evolution and the future singularities for this 
model. Compared with the standard model of cosmol- 
ogy, this model has had three additional parameters, Co, 
Ci and po : choices of £o and negative po can cause the 
universe accelerate; (1 can drive the cosmic fluid from 
the quintessence region to the phantom one [5j , and pos- 
itive po may both prevent the big rip for phantom phase 
and lead to the Type III singularity of Ref. |(| for the 
quintessence phase. The relation between the choices of 
parameters and the future singularities of the cosmologi- 
cal evolution in this extended model is summarized as in 
the following table and we expect more detail investiga- 
tions on viscosity effects to be carried out. 
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Parameters 


Future singularity 








(at t — ► t s ) 


7 < 


Po < 


a — > oo, p — > oo 




Po > 


1 + 7 6> T - T/Ti > 


No 






1 + 7 6> T - T/Ti < 


a — ► oo, p — ► oo 


7 = 


No 


7 > 


po < 


No 




Po > 


1 + j6 T - T/T x > 


No 






1 + 7 O T - T/Ti < 


a — > 0, p — * oo 
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